Fidelity recovery in chaotic systems and the Debye- Waller factor 
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Using supersymmetry calculations and random matrix simulations, we studied the decay of the 
average of the fidelity amplitude / £ (t) = (V'(O)I exp(27ri-ff e T) exp(— 2mHoT)\ip(0)), where H t differs 
from Ho by a slight perturbation characterized by the parameter e. For strong perturbations a 
recovery of / e (r) at the Heisenberg time r = 1 is found, ft is most pronounced for the Gaussian 
symplectic ensemble, and least for the Gaussian orthogonal one. Using Dyson's Brownian motion 
model for an eigenvalue crystal, the recovery is interpreted in terms of a spectral analogue of the 
Debye- Waller factor known from solid state physics, describing the decrease of X-ray and neutron 
diffraction peaks with temperature due to lattice vibrations. 
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The concept of fidelity has been developed by Peres 
as a tool to characterize the stability of a quantum- 
mechanical system against perturbations 0. It was in- 
troduced as the squared modulus of the overlap integral 
of a wave packet with itself after developing forth and 
back under the influence of two slightly different Hamil- 
tonians. Very similar concepts had been applied already 
in the old spin-echo experiments of nuclear magnetic res- 
onance half a century ago (see reference for a review). 
The renewed interest in the topic results from the idea 
to realize quantum computers by means of spin systems, 
where stability against quantum-mechanical perturba- 
tions obviously is of vital importance [2| . 

Roughly speaking there are three regimes. In the per- 
turbative regime, where the strength of the perturbation 
is small compared to the mean level spacing, the decay 
of the fidelity is Gaussian. As soon as the perturbation 
strength becomes of the order of the mean level spacing, 
exponential decay starts to dominate, with a decay con- 
stant obtained from Fermi's golden rule 0, . For very 
strong perturbations the decay becomes independent of 
the strength of the perturbation. Here, in the Loschmidt 
regime, the decay is still exponential, but now the decay 
constant is given by the classical Lyapunov exponent [£j . 
Exactly such a behaviour had been observed experimen- 
tally in a spin-echo experiment on isolated spins coupled 
weakly to a bath of surrounding spins Q • 

Gorin et al. || calculated the fidelity decay within 
random matrix theory in the regime of small perturba- 
tions and could correctly describe the change from Gaus- 
sian to exponential behaviour with increasing perturba- 
tion strength. The Lyapunov regime is non-universal and 
thus not accessible in a random matrix model. 

Intuitively, one would expect that in chaotic systems 
the fidelity decay is stronger than in integrable ones. The 
opposite is true. Prosen et al. showed that a chaotic 
system is much more fidele than a regular one, and sug- 
gested to use chaotic systems in quantum computing to 
suppress chaos. 

It will be shown here that this is only part of the truth, 



and that for chaotic systems there is even a partial re- 
covery of the fidelity at the Heisenberg time. This work 
extends the results by Gorin et al. [S| to the regime of 
strong perturbations using supersymmetry techniques. It 
is stressed that our results are generic and not restricted 
to random matrix systems. In a spin-chain model, e.g., 
the fidelity recovery has been observed recently as well 
0.. 

Using the Brownian-motion model for the eigenvalues 
of random matrices introduced by Dyson many years ago 
it will be shown that this behaviour has its direct 
analogue in the Debye- Waller factor of solid state physics. 
We shall sketch the calculation for the GUE only, and 
will just cite the result for the GOE. More details will be 
presented in a forthcoming paper 

Let us start with an unperturbed Hamiltonian Hq, 
which is turned by a small perturbation into 



i?0 = Hq cos 4> + Hi sin < 



(1) 



It is assumed that both Hq and H± have a mean level 
spacing of one, i. e. the variance of the off-diagonal ele- 
ments is given by N/ ir 2 0] . This particular choice of the 
perturbation guarantees that the mean density of states 
does not change with e. It will be come clear below that 
(j> has to be of 0(1/ VN) to allow a well-defined limit 
N — > co. We therefore introduce e = 4N(f) 2 as the per- 
turbation parameter. It follows for the fidelity amplitude 



Mr) = (V(Q) 



(2) 



where ip(0) is the wave function at the beginning, and 
c = cos <p, s = sm(f). The squared modulus of / c (t) yields 
the fidelity F e (r), originally introduced by Peres The 
calculation of the average of F £ (r), however, is technically 
more involved and therefore not considered here. 

For the present work the paper of Gorin et al |8| is 
of particular relevance. The author considered a slightly 
different parameter dependence 



H x = H + AV , 



(3) 



2 



where the variance of the off-diagonal elements of V was 
assumed to be one. Comparison with equation JQ| shows 
that the respective perturbation parameters are related 
via A = y/e/ (2tt). 

In the paper by Gorin et al. [8j the Gaussian average 
of the fidelity amplitude was calculated in the regime of 
small perturbation strengths, correct up to 0(e), 



</e(r))~e 
where C(r) is given by 



-eG(r) 



(4) 



C{t) = j + T -- £ J\(t')dt'dt, (5) 

and b 2 (r) is the spectral form factor. f3 is Dyson's univer- 
sality index, i. e. (3 = 1 for the Gaussian orthogonal en- 
semble (GOE), (3 = 2 for the Gaussian unitary ensemble 
(GUE), and (3 = 4 for the Gaussian symplectic ensemble 
(GSE). Equation J3J describes correctly the change from 
Gaussian to exponential decay with increasing perturba- 
tion strength. Equation |@J was the motivation for the 
introduction of e as the perturbation parameter. 

In chaotic systems the average of the fidelity amplitude 
over the initial state ^(0) reduces to a trace, 



(/e(r)) = ^(Tr 



(/ c (r)) may be written as a Fourier transform, 
(/ e (r)) = J dE 1 dE 2 e 2 ^- E ^R t (E 1 ,E 2 ) 



(6) 



(7) 



where 



Re (El,E 2 ) 

1 



N 



Tr 



Ei- - cHo - sHx E- 



2+ 



Ho 



(8) 



with E± = E ± i?/. Using standard supersymmetry tech- 
niques |14|. this can be written as 

Re (El,E 2 ) ~ 

— J d[x] d[y] J2( X *n x m - ^m){y* m Vn ~ V*mVn) 



xe 



x ( e 



-{[x+Bix-ytBay] 
^[cx+ffox-ytffoy] 



(9) 



where x = (a^fi . ..x N ,£, N ) T , y = (yi,m . ..y N ,VN) T , 
and 

d[x] = Y[ dx n dx* n d£ n dC n , d[y] = Y[ dy n dy* n d-q n drf n . 

n n 

(10) 

We adopt the usual convention using latin letters for 
commuting, and greek ones for anticommuting variables, 



respectively. Now the Gaussian average over the matrix 
elements of Hq and H\ can be performed elementary. 

The subsequent steps (Hubbard-Stratonovich transfor- 
mation, integration over the x, y variables, saddle point 
integration etc.) are essentially the same ones as for the 
calculation of the spectral form factor (see e. g. chapter 
10 of referen ce I l5il for the GUE case). Details can be 
found in Ref. [l2l We then arrive at the result 

oo 1 

d it? j? \ p 2 f a I a x + y + x 2 -y 2 

Re{Ei,E 2 ) ~ -JdxJ dy {x + y)2 



y , e ~2TTipE(x+y) e -^p 2 (x+y)(l+x-y) |-qs 



where p = y/l - [ttE/{2N)] 2 is the mean density of 
states, and E = {E\ + E 2 )/2, E = E\ — E 2 . Inserting 
this result into equation (JJJ), and introducing E and E 
as new integration variables, we get, fixing the constant 
of proportionality by the condition / c (0) = 1, 

oo 1 


xS[t - p(x + y)] e --2p 2 ^+v)^-y) . (12) 

The E integration is nothing but an energy average. In 
the limit TV — > oo only the band centre contributes where 
p takes the constant value one. We may then discard this 
average and obtain 



(/e(r)> 



1 



Min(T,l) 



dy(l + T -2y)e 2 



r(l+r-2y) 



(13) 



The integral is easily performed with the result 
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where 
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s{x) 



sinh(x) 



T < 1 

r>l ' 
(14) 

(15) 



We have thus obtained an analytic expression for the 
GUE average of the fidelity amplitude for arbitrary per- 
turbation strengths. In the limit of small perturbations 
it is in complete accordance with the result obtained by 
Gorin et al. 

The calculation for the GOE is done in exactly the 
same way. It is technically much more involved, but for- 
tunately most of the work for this case has already been 
done by Verbaarschot, Weidcnmullcr, and Zirnbauer in 
their disseminating work |14| . In this way we get for the 
GOE average of the fidelity amplitude 



(/e(r)> =2 J du 

Max(O.T-l) 



v dv 



y/[u 2 ~V 2 ][(u+l) 2 -V 2 ] 



3 



(t-u)(1-t + u) 

(V 2 - T 2 ) 2 

x[(2m+1)t- 



(16) 



T 2 + w 2 l e -54(2«+l)r-r 2 +t , 2 ] 



To affirm the analytical findings, random matrix sim- 
ulations were performed for all Gaussian ensembles in- 
cluding the symplectic one, which has not been treated 
analytically. The average was taken over up to 8000 ran- 
dom matrices of rank TV = 500 for H , and for each of 
them over 50 random matrices for Hi. 

The results are shown in Figure ^ For the GOE and 
the GUE the numerical simulations are in perfect agree- 
ment with the analytical results for all e values shown. 
For comparison, the fidelity amplitudes in the linear re- 
sponse approximation J3J are presented as well. For small 
perturbation strengths or small values of r the linear re- 
sponse result is a good approximation, but the limits of 
its validity are also clearly illustrated. 

For small perturbation strengths e the decay of the fi- 
delity is predominantly Gaussian which changes into a 
behaviour showing a cross-over from an exponential to a 
Gaussian decay at e w 1, in accordance with literature. 
The most conspicuous result of the present letter, how- 
ever, is the partial recovery of the fidelity at the Heisen- 
berg time r = 1 which has not been reported previously 
to the best of our knowledge. This recovery is most pro- 
nounced for the GSE and least for the GOE. This is il- 
lustrated in Figure |21 showing a direct comparison of the 
fidelity amplitudes of the Gaussian ensembles for e = 30. 
Preliminary simulations showed that the fidelity recovery 
is observed not only in the fidelity amplitude, but in the 
fidelity as well. 

What is the origin of the surprising recovery? We 
believe that there is a simple intuitive explanation in 
terms of Dyson's Brownian motion model |ll| . Since 
the mean density of states is kept constant during the 
parameter change, the eigenvalues of H^ may be written 
as = k + 8%, where 8% fluctuates about zero. For 
strong perturbations the eigenvectors of the perturbed 
and unperturbed system are uncorrelated, and we get 



</^)>^E<i^i> 27riT(fc ^ 



(17) 



A7 



where the Rki are obtained elementary from the eigen- 
vectors, and W is given by 



W 



-(2-r) 2 (A- 2 > 



(18) 



In the second step a Gaussian approximation was ap- 
plied. Within the framework of the Brownian motion 
model (<5 2 ) is interpreted as the mean squared displace- 
ment of an eigenvalue from its equilibrium position. It 
is proportional to "temperature" T, which is just the re- 
ciprocal universality factor /3, whence follows 
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FIG. 1: Average of the fidelity amplitude (/ e (r)> for the GOE 
(a), the GUE (b), and the GSE (c) for e = 0.2, 1, 2, 4 and 
10. t is given in units of the Heisenberg time. For the GOE 
and the GUE the solid lines show the results of the analytical 
calculation, for the GSE of the numerical simulation. The 
dotted lines correspond to the linear response approximation 
0. The numerical results are in agreement with the analyt- 
ical results within the limits of the line strength. 



with some constant a. It follows from equation (|17fl that 
there is a revival of the fidelity at the Heisenberg time r = 
1 decreasing with "temperature" proportional to e~~ aT . 
This is exactly the behaviour illustrated in Figure [3 

There is a perfect analogy to the temperature depen- 
dence of X-ray and neutron diffraction patterns in solid 
state physics. Caused by lattice vibrations the intensi- 
ties of the diffraction maxima decrease with temperature 



4 




FIG. 2: Average of the fidelity amplitude (/e(t~)) for e — 
30. The dashed and dashed-dotted lines show the analytical 
results for the GOE and the GUE, respectively. The solid 
fine corresponds to the numerical simulations for the GSE, 
reliable up to 10 -3 . 



with a dependence described by the Debye- Waller factor 
W DW = e- pg2T , (20) 

where (3 is another constant, and g is the modulus of the 
reciprocal lattice vector characterizing the reflex (see e. g. 
appendix A of reference This is our justification to 

call W a spectral Debye- Waller factor. 

One may argue that due to equation (|17fl there should 
be revivals for all integer multiples of the Heisenberg 
time, which are not observed. This can be understood by 
considering the analogy between the spectral form factor 



(21) 



and the structure factor in condensed matter, 

S (k) = I ^ e 2^.(fl n -J? m ) (22) 



where the R n are the positions of the atoms, and k is 
a point in the reciprocal lattice. In liquids and glasses 
these peaks are smeared out, and the structure factor- 
depends only on k = \k\. S(k) starts at zero for k = 0, 
climbs up to two to three at the k value corresponding 
to the reciprocal atomic distance, oscillates about one 
with decreasing amplitude for lar ger k values, and ap- 
proaches one in the limit k — > oo |17| . There is a strik- 
ing similarity with the spectral form factor K{t) of the 
GSE, the eigenvalue "crystal" with the lowest available 
temperature. K(r), too, starts at zero for r = 0, has 
a logarithmic singularity at the Heisenberg time t = 1, 
and approaches one for r — * oo. The oscillations found 
in S(k), however, are absent in K (r). Fidelity amplitude 



and spectral form factor thus do not show any structure 
at multiple integers of the Heisenberg time, in contrast 
to the corresponding condensed-matter quantities. The 
explanation is straightforward: Pauli's principle prevents 
the atoms from approaching too closely, whereas in the 
eigenvalue "crystal" such a lower limit does not exit. 

If it is true that the universal spectral properties for all 
chaotic systems are described correctly by random ma- 
trix theory [isj . and there is overwhelming evidence for 
this fact, then the exact expressions for the fidelity am- 
plitude decay derived in this paper are universal and hold 
for all chaotic systems. But this means that the fidelity 
decay in this regime reflects the phase space properties 
of the system, but not its stability against perturbations 
as originally claimed by Peres Q. Only in the regime of 
the Loschmidt echo the system stability shows up. 

This paper has profited a lot from numerous discus- 
sions on the subject of fidelity with Thomas Seligman, 
Cuernavaca, Mexico, Thomas Gorin, Freiburg, Germany 
and Tomaz Prosen, Ljubljana, Slovenia. Thomas Guhr, 
Lund, Sweden is thanked for discussions of the supersym- 
metry aspects of this paper. The work was supported by 
the Deutsche Forschungsgemeinschaft. 
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